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Abstract 

The energy in a square membrane Q subject to constant viscous damping on a subset 
u G Q decays exponentially in time as soon as u satisfies a geometrical condition known 
as the "Bardos-Lebeau-Rauch" condition. The rate t(lo) of this decay satisfies t{uj) = 
2min{—iJ.{ui), g{uj)) (see Lebeau [9]). Here /i(tt') denotes the spectral abscissa of the damped 
wave equation operator and g{u!) is a number called the geometrical quantity of u and defined 
as follows. A ray in is the trajectory generated by the free motion of a mass-point in Q 
subject to elastic refiections on the boundary. These refiections obey the law of geometrical 
optics. The geometrical quantity is then defined as the upper limit (large time asymp- 
totics) of the average trajectory length. We give here an algorithm to compute explicitly g{uj) 
when w is a finite union of squares. 

MSG Numbers: 35L05, 93D15 
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Let = [0, 1]^ be the unit square of and let C 17 be a subdomain of ^. We are interested 
here in the problem of uniform stabilization of solutions of the following equation 

o 

uuix,t) - Au{x,t) + 2xui{x)ut{x,t) = 0, a; £17, i > 0, 

u{x,t) = 0, a; e do., t>0 , . 

uti.,0)^ui e L2(r!), 

where x^j is the characteristical function of the subset uj. This equation has its origin in a physical 
problem. Consider a square membrane fi. We study here the behaviour of a wave in fl. Let u{x, t) 
be the vertical position of x G 17 at time t > Q. We assume that we apply on w a force proportional 
to the speed of the membrane at x. Then, u satisfies equation (jO.ip . To get more information on 
this subject, one can refer to [7], [TU], With regard to the one-dimensional case, the reader can 
consult [S] and [3]. 
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We do not deal here with the existence of solutions. We assume that there exists a solution u. Let 
us define the energy of u at time t by 

E{t) = / \\Iu\^ + uldx. 
Jn 

It is well known that, for any t > 0, 

E{t) < CE{0)e-^^\ (0.2) 

where C, t > 0. As proven by E. Zuazua [13], this result remains true in the semilinear case when 
the dissipation is in a neighborhood of a subset of the boundary satisfying the multiplier condition. 
Let us now define 

Definition 0.1 The exponential rate of decay t{uj) is defined by 

t{uj) = sup{t > s.t. 3C >0 for which fO|) holds }. 

Many articles have been devoted to finding bounds for r(w). The reader can consult [T], [5] and 
[llj . In the case of a non-constant damping, the reader may see j4]. G. Lebeau proved in [9] that 

t{lu) = 2 min(— (7(w)). 

Here fJ.{io) denotes the spectral abscissa of the damped wave equation operator and g{u) is a 
number called the geometrical quantity of w and defined as follows. A ray in Q is the trajectory 
generated by the free motion of a mass-point in ^l subject to elastic reflections in the boundary. 
These reflections obey the law of geometrical optics: the angle of incidence equals the angle of 
reflection. If a ray meets a corner of the reflection will be the limit of the reflection of the rays 
which go to this corner. It is easy to verify that the ray runs along the same trajectory before and 
after the reflection but in an opposite direction. 

Let 7 : ^ £7 be the parametrization of the ray by arclength. Let 

C = { rays in fi} . 

Furthermore, let po = (Xq, a), where Xq € and a G [0, 2n[. Consider the ray which starts at 
Xq in the direction of the vector (cos(a), sin(a)). Each ray of C can be deflned in this way. More 
precisely, if F = x [0, 27r[ then 

C = {ipo s.t. Po e F}. 

In the whole paper 7 G C will be noted 7 = [Xq, a]. If 7 G C and if t > is a positive real number, 
we write: 

It = l/[o,t]- 

Let 

t>0 

be the set of paths of finite length not necessarily closed. For jt € C of length < > 0, we define: 

1 /■* 

"^(7t) = 7 / Xu,{'yt{s))ds. 
t Jo 
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The definition of m can be extended to C. Indeed, if 7 S C, set 

to(7) = limsupm(7t). 

t — *+oo 

It is proven in Section 1.2 that 

771(7) = li™ '^{it)- (0-3) 

t — y-\-oo 

For a ray 7 belonging to C, 771(7) represents the average time that 7 spends in 17. 
Definition 0.2 The geometrical quantity g{uj) is defined by 

g{uj) = Km sup inf 777(7^). 

t— >+cxD tsC 

We are interested here in a precise study of the geometrical quantity g. The first part of this paper 
is devoted to studying the rays. We recall some well known properties of rays. In the second part, 
another expression for g is given. Namely, we prove that ii lj G ^ whose boundary is a finite union 
of C^-curves, then 

g{io) = inf limsup777(7f) = inf 777(7). 

It is easier to work with this second definition of g{Lij). An important application of this theorem 
is given in the third part: we obtain an algorithm which gives an exact computation of g{u}) when 
w is a finite union of squares. This work has many interests. At first. Theorem 12.11 gives another 
definition for 3(0;), much more easy to manipulate than the original one. Secondly, maximizing 
the exponential rate of decay is interesting from the point of view of physics. For these questions, 
knowing exactly g{ijj) is important. We give some exact computations of 3(0;) with the help of our 
algorithm at the end of the paper. An interesting question is then: how can we choose the subset 
uj such that g{Lj) is maximum? At the moment, this problem is still open. The following inequality 
is always true g{ijj) < 1^1 (see Section 2) where |a;| is the area of oj. Let us set for a G [0, 1] 

S(a) — sup giiu). 

\Lj\—a 

Some questions then arise naturally 

- What is the value of 5(a)? In particular, is it true that S{a) = a? 

- Can we find a subset w for which g{uj) — \lli\7 

- Can we find an optimal uj (i.e. an lo that maximizes g) among the domains that satisfy the 
multiplier condition? 

If \uj\ = or 1^1 = 1, the answers are obvious. However, if \uj\ e]0, 1[, these questions seem to be 
much more difficult and are still open. 

1 Basic properties of rays 

1.1 Different representations of rays 

In this section, we regard billards in from different point of views (see for example [H]). In- 
stead of reflecting the trajectory with respect to a side of fi, one can reflect with respect to 
this side. We then obtain a square grid and the initial trajectory is straightened to a line. This 
gives a correspondance between billard trajectories in and straight lines in the plane equipped 
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with a square grid. Two lines in the plane correspond to the same billard trajectory if they differ 
by a translation through a vector of the lattice 2Z + 2Z. The factor 2 in 2Z + 2Z is important. 
Indeed, two adjacent squares have an opposite orientation in the sense that they are symmetric 
with respect to their common side. 




1 2 

Figure 1: Point of view of infinite plane 

Now, consider T = [0, 2f and identify its opposite sides. A trajectory then becomes a geodesic 
line in the flat torus. 




1 2 



Figure 2: Point of view of flat torus 

Finally, trajectories in O can be seen in three different ways 

- the definition which consists in reflecting trajectories in sides of fl 

- the point of view of straight lines in 
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- the point of view of geodesic lines in a flat torus T. 

Note that the set w must be reflected in the same way as we did for Vt. This means that the sets 
(jj contained in two adjacent squares must be symmetric with respect to their common side (see 
Figures [Hand HI). 

Let t> Q and 7 G C. The value 

XMs))ds 

does not depend on the point of view we adopt. 



1.2 Open and closed rays 

A ray is said to be closed if it is periodic in the flat torus. It is said to be open if it is not closed. 
We have the following characterization. The ray [Xq = (xq, yo); ct] is open if and only if cos a and 
sin a are independent over Z, i.e. if and only if tana S K — Q. Indeed, for example using the point 
of view of a flat torus, the trajectory is periodic if and only if there exist < ti < t2 and rt, m G N 
such that 

J xo + ti cos a = xo + ^2 cos a-\-2n 
\ 2/0 + ^1 sin a = yo + ^2 sin a + 2m 

Hence 2m{ti — i2)cosa — 2n{ti — t2)sma = 0. Let us now consider open rays. They have the 
following properties (see for example [5] p. 172). 

Proposition 1.1 Let 7 = [Xo,Q!] be an open ray. Then, {74, t > 0} is dense in the torus T and 
then in fJ. In addition, if uj is quarrable (i.e. Xui Riemann-integrable onfl) then m{pit) > l^^^l- 

t — *+oo 

More precisely, let a be such that tana G K — Q then 

Ve > 0, 3to 1 1 > to, VXo G VL : ifjp^ = [Xq, a] then \m{jt) - |w|| < e. 
In other words, t^ is independent of Xq . 

Here, |w| stands for the area of Consider now closed rays. They satisfy the following properties: 

Proposition 1.2 Let 7 = [Xo,a] be a closed ray. Then, there exist p,q G Z relatively prime such 
that: 

P , . q 

cos a = — j^=^= and sm a = — j^=^= . 

vp^ + q^ Vp^ + q^ 

The period L of the trajectory is L = + . We have 

Remark 1.1 The assertion \0.3\) follows immediately from Propositions and^TT^ Moreover, 
let "f £ C. If "f is open, we have 771(7) — 1^1 '^'^'^ '/ 7 periodic of period L, we have 7*7,(7) = 
x/o xU^{s))ds. 

As one can check, an immediate consequence of Proposition 11.21 is the following result 
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Proposition 1.3 Let uj C Q be Riemann-integrable, then 

g{^) < 1^1- 

Closed rays can easily be described. Let us adopt the point of view of a flat torus. Let us consider 
a ray 7 — [Xo,a] with Xq E T and cos a — p/ \/jP + q^, sin a = q/\/p'^ + q^, P and q relatively 
prime. We assume that p>0 and q>0 (i.e a £ [0, 7r/2]). From remark [5TT] below, the study can 
be restricted to such rays. 

Theorem 1.1 The ray 7 is an union of p + q parallel segments directed by u — cosai + sinaj. 
Among these segments, p of them start at Ak, 1 < k < p, and the q others start at B^, 1 < k < q 
with 

Ak (^0, {yo + q/p{2k - xo) - 2E{^{yo + q/p{2k ~ xo)))^ 

and 

Bk {{x^ + p/q{2k ~ yo) - 2E{^{xo + p/q{2k - yo))) , 0^ ). 

In addition, the distance between two neighbouring segments is constant and equal to 6 = 4/L, 
where L = 2^ p^ + cp- is the period. 
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Figure 3: Closed ray for p = 3 and g = 2 in the infinite plane 

Proof : In the infinite plane R^, the ray is defined by x(t) — xq + 2tp/L and y{t) = yo + 2tq/L, 
where L — 2^p^ + q^ is the period. Let us restrict the study to the segment [XoXi] of length L. 
In ]XoXi[, the ray intersects p vertical straight lines defined by the equations x = 2k, fc = I, • • • ,p 
and q horizontal straight lines defined by the equation y ~ 2k, k — 1, - ■ ■ ,q. This gives p + q 
segments starting at points with coordinates 

(0, yo + <7/p(2A: - Xo)) 

and 

(a;o+p/g(2fc-yo),0). 
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Figure 4: Closed ray p = 3 and q = 2 in the flat torus 



Coming back to the point of view of the torus, this gives p + q segments starting at Ak and Bk ■ 
Since p and q are relatively prime, it is well known that {kq \p], 1 < k < p} = {0,1, . . . ,p — 1}. 
Hence, if A^^ and Ak2 are two neighbouring points, Aj-^Aj-^ = 2/p. As a consequence, the distance 
between the straight lines ( , u) and {A^^ , u) is 



AkiAi~2-u 



2 

-J 



-qi+pj 



P + q2 L 



In the same way, if Bk^ and Bk2 are two neighbouring points Bk^Bk2 =2/q, the distance between 

{Bki,u) and {Bk^^u) is 

^ ^ 4 
1' 



Bk^Bk^.u 



2^_qi--pj_ 
9 ' + 9^ 



Let fci, ^2 be such that 



dist{Aki,0) = mmdist{Ak,0) and dist{Bk2,0) = m.mdist{Bk,0). 

k k 

It remains to prove that the distance between the two segments starting at Ak^ and Bk2 is 6 = 4/L. 
Since the distance between two points Ak is 2/p, there exists n e N such that the ordinate y* of 
Aki satisfies 

2n<ya + q/p{2ki - Xq) - 2E{^{yf, + q/p{2ki - x^)) <2n + 2/p 

=^ y* =yo + q/p{2ki - Xq) - 2E{^{yo + q/p{2ki - Xq)) - 2n < 2/p. 
In the same way, there exists m G N such that the abscissa x* of Bk2 satisfies 
2m < Xo+p/q{2k2 - yo) - 2E{^{xo +p/q{2k2 - yo)) < 2m + 2/q 

=^ X* = Xo+p/q{2k2 - yo) - 2E{^{xo+p/q{2k2 - yo)) - 2m < 2/q. 
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The distance S between the two segments verifies: 

\qx* +py* 



S 



Ak^Bk^.u-^ 



An easy computation shows that \qx* + py*\ is an even integer number. Consequently, S can be 
written as Ak/L, where k is an integer. Since < x* < 2/q and < y* < 2/p, we get 6 < 8/L, 
and hence, S = A/L. □ 

Remark 1.2 // a quarrahle domain u is such that g{uj) — \uj\ then for all a G [0,7r/2] and for 
almost all Xq G T, m(7) = \uj\, where 7 = [Xo,a]. 

Indeed, consider an open ray 7. Then 771(7) = l"^!- Let Xq G Assume that a ray 7 starting at 

> 

Xq is closed of length L. Let M and N be two points of 7 such that AIN is orthogonal to U and 
MN = 6. Let also P{s) = sM + (1 - s)N and 7^ = [^'(s), a], then 

5r(tj) < inf to(7£) < / m{-fl)ds < \uj\. 

If .g(ti-') = then m{jf^) — \uj\ almost everywhere. If duj is a finite union of curves of class C^, 
then the number of discontinuity points of m(7£) is finite. It is null if duj does not possess any 
segments. The remark above immediately follows. 



2 An equivalent definition for g 

2.1 The geometrical quantity g' 

We define, for any set ui: 

Definition 2.1 The geometrical quantity g' {uj) is defined by 

g'{uj) = inf to(7). 

As easily seen, g' is easier to study than g. As shown in Section 3, g'iio) can be computed explicitly 
for a certain class of domains oj. Together with Theorem 1 2. 1[ this gives an algorithm for computing 

Theorem 2.1 Let uj be a closed set whose boundary is a finite union of curves of class . Then 

g{uj) = g'{u}). 

As a first remark, the same argument as in the proof of Proposition 11.31 shows that 



0<g'H<k|. (2.1) 

The proof of Theorem 12. II is given in the appendix. It is easy to see that (7(1^) < g'{uj)- Inequality 
g'i'-^) ^ 3(^) is much more difficult to obtain. To prove Theorem 1 2. 1[ we consider a sequence of 
rays 7" — [a;„,0„] and a sequence of real numbers tn — > +00 for which 771(7^") — > g{uj). After 
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choosing a subsequence, there exists {x,9) e ft x [0, 27r[ such that hm„ a:„ = x and hm„ = 6. 
We will now show that lim„ m(7^) — m{[x, 0]t^) — 0. The conclusion then follows. The difficulty 
in this proof is that the function {x,6) m{[x,6]t) is not continuous. A direct consequence of 
Theorem l2.1l is that it suffices to consider closed rays in the explicit computation of g{uj). Namely, 
let Cc be the set of closed rays, then we have: 



g{u;) = inf m(7). 



(2.2) 



3 Explicit computation of g 

We prove in this section that for a particular class of domains lu, the properties of the geometrical 
quantity g we obtain above allow to compute explicitly g{Lu). Let N € N* , and let a; C 51 be a 
finite union of squares {Cij)i<ij<N , with 



Ci. 



1 i 



N N 



J - f J_ 

N ' N 



Obviously, u is Riemann-measurable, closed and its boundary is a finite union of curves. 



3.1 Influence of a 

The first result we obtain is the following 

Theorem 3.1 Let 7 = [Xq, a] be a closed ray such that tana = q/p, with p and q relatively prime 
and p + q > 27V, then 

iV2 

I m{j) - |w| I < mm{|w|, 1 - \uj\} < —. 

pq 2pq 

Proof : We adopt the point of view of flat torus T = [0, 2]^. At first, let K be one of the Cij i.e. 



K = a, 



^0 - 1 io_ 
N ' iV 



Jo - 1 Jo 
N ' N 



ABCD. 



The period of 7 is i = 2y/jp + q^. Let S ~ A/L. Then cos a = p6/2 and sin a = qS/2. By 
Theorem II. 1[ the ray 7 is constituted by p + q segments directed by u. The distance between two 
neighbouring segments is 5. 



Let us define the orthonormal frame TZu = (X, /„, J^), where X is such that XD.u — 0, XD.XA ■ 
0, = Yd/\\XD\\, and X = 'XA/\\YA\\ = u (see Figure^. 



In this frame, the coordinates of A,B,C and D satisfy 

, sin a + cos a cos a . 



. sma, „,smQ! sma + cosa. 



N 



N ' 



iV 



C{ 



N 



N 



The ray 7 is constituted of p + g vertical segments. R of them meet K (see Figure |6]). We note 
a;o,xi, . . -XR-i their abscissa 

a;o = T' e [0, 8\ 

Xi — xq + iS, < i < R — I 
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1 



2 



Figure 5: Ray 7 = [0, a] with tana = 4/15 



Since a+cos a _ ^ there exists at least one segment in K and R>1. 
Let 



and 



Then 



and 



/2W 



sm a sm a sm a + cos a , , 
cot a x-\ = cot a [x — — ) H — line (ABj 



N 



N 



N 

, sm a , sm a + cos a 
■tana(a;-^) + 

, sin a + cos a, cos a 
■ tana \x — ) ^ — Ime \BC) 



N 



N 



, . sma , cos a. . , 
gi{x) =— tanaxH ;^ = — tana(a; Ime (AL>) 



52(2;) 



N ' N 
. cos a , , sm a + cos a , cos a / ^ „n 
cot a (a; — — j=cota(a; — )H — — Ime (CL'). 



N 



N 



N 



N 



if a; < 



fi{x) ifO<x< 



sm a 

/2(x) if ^ < x < ^ ^ 

cos a -f ^ sin g+cos a 
Jy 



N 



if X < 



cos a 

N 

gi{x) ifO<x< 

N 



cosg 



52(2;) if Sep < a; < "'""+^°^" 



cos a 



rQ Sin g+CQS a 



N 



Jr Jo 



{f{x)-g{x))dx = \K\ 



1 



1 s 



Let k be the number of vertical segments which meet [AB] , i.e. the k segments such that 



r + {k - 1)6 < — ^ < r + fcJ. 
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Figure 6: Vertical segments in the square ABCD 



If 1 < fc < i? — 1, if at least one segment meets [AB] and if at least one segment meets [BC] then, 
as one can see 



i—O i—0 i—k 



■.,+S/2 



t=0 •'Xi-S/2 
^ ,r+(k-l/2)5 

4 Jr-S/2 



fi{x)dx + - 



k -x,-&/2 
^ l.r+{B-l/2)S 



r+(fe-l/2)<5 



f2{x)dx 
f2{x)dx. 



If k ~ or k ~ R, equation (|3.ip remains true. Hence, in all cases 

R-l 

f(xA ^ - " ' f(x)dx - ^ 



(p+q)S/2N ^ |.r-A72 

f{x)dx - - fi{x)dx 
4 Jo 



1 <.r+(fe-l/2)<5 -, i-r+(R-l/2)5 

{h{x) ~ h{x))dx + - h{x)dx. 



qS/2N 



{p+q}S/2N 



In the same way, let I be the number of segments which meet [AD] 

cos rv 

r + {l-l)6< -j^ <r + l6. 



(3.1) 
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We obtain that, in all cases 



R-l 



^ rT+{l-l/2)S 



gi{x)dx + 



^ ,-r+{R-l/2)5 



i=0 



r-5/2 4 Jr+{l-l/2)5 

(p+q)6/2N nr-il2 

g{x)dx - -^j gi{x)dx 



g2{x)dx 



+ 



1 



r+((-l/2)5 ^ 

{9i{x) - g2{x))dx + - 

pS/2N * 



r + (/f-l/2)(5 



g2{x)dx. 



{p+q)S/2N 



Finally 



4Ar2 



r-S/2 . 

h-91 + 7 



1 

'i 
+ 

'pS/2N 

tan a + cot a 



r+(k-l/2)S 



/1-/2 



qS/2N 



^ /.r+(i-l/2)5 



r+(fi-l/2)(5 

91-92 + -; I f2-92 



4/, 



4 J{p+q)5/2N 



[-(r - 6/2f + {r + {k- 1/2)6 - q5/2Nf 

o 

+(r + - 1/2),^ - p5/2Nf -{r+{R- 1/2)5 - (p + q)S/2N)'^] . 

However, in this last equality, each of the four square terms is less than 6/2 and tana + cot a = 
1 4 L2 

= = hence 

sin a cos a pq6'^ Apq 



XK{l{s))ds 



4Af2 



< 



4pq' 



(3.2) 



We now consider: 



where £^ is a subset of {1,2,..., 2Ny and card(£;) = N^\u}t\ = 4:N^\u}\ ; then 



{i,j)eE 



Hence 



This leads to 



Ml) -H= ( 7 / xc.,. 



{l{s))ds 



4Ar2 



, , , , , . ^ card(£;) , 
I m(7) - I < — . = — \uj\ 



4pq pq 

If the area of oj is greater than 1/2, let oj' be the adherence of O — a; which is also a finite union 
of squares Cij , Hence 



X^'{l{s))ds- il-\uj\) 



<—{i-H 

pq 
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or L L 

Finally, we get 1 772(7) - I < (1 - H)- □ 

pq 



3.2 Representation of u 

For convenience, consider the frame (O, /, J) with / = i/{2N) and J = j/ {2N). As a consequence, 
il is now the square [0, N]'^, and T is the torus [0, 27V] ^. The length of closed rays defined by (p, q) 
{p,q relatively prime) such that cos a = p/^/p'^ + q^ and sin a = qj \/p^ + q^ must be modified. 
The period is now 

L = 2iVv/p2 + g2 

and Ci J is such that: 

C»j = [« - X [j - 1, j]. 

Let us define the matrix M N x N such that Af^j = 1 if Cj^i C c^j, and Mij = if Cj,i {Z! With 
this notation, the matrix M is a representation of cj. As an example, if oj is as in Figure [7] then 



M = 



/ 1 


1 








1 


1 \ 


1 


1 


1 


1 





1 





1 





1 





1 














1 


1 


1 


1 














I 


1 


1 








/ 



"1 



1 2 3 4 5 6 
C C C C C C 

1 2 3 4 5 6 

columns 



Figure 7: A domain uj for iV = 6 
Let us also define the 2A^ x 2N matrix A such that Aij = 1 if Cj,i C ojt, = if Cj^i <^ lut, i-e. 



A. 



1,3 



A 



2N+l-i,j 



A 



i,2N+l-j 



A 



2N+l-i,2N+l~j ■ 



(3.3) 
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3.3 Influent rays 

Let 7 be a horizontal ray starting at Xq — (xqjI/q). Then, 

N 



Assume that ya — k is an integer. Since each Cij is closed, 

N 

N 



1 ^ 



Obviously, this ray does not realize the infimum in the definition of g. In the same way, for a 
vertical ray, only the numbers 

1 ^ 

— VM„-, l<i<N. 

must be considered. Now, let us deal with oblique rays and fix a. Let us denote tana by g/p, 
where p and q are relatively prime integers. 

Remark 3.1 One has to consider only the rays defined by ^ = [Xo,a\ with Xq G T and a G 
[0,7r/2], i.e. p>0 and q>0. 



2 



1 2 

Figure 8: Starting points with a G [0,7r/2] 

Indeed, keeping the notations of FigurelH the ray 7^ — [Xo,ui] is equivalent to 7^ — [Xi,u\. In the 
same way, the rays 7^ — [Xq, U2] and 73 = [Xq, U3] are equivalent to 7^ = [-''^2, w] and 7'^ — [X3, u\. 

Let us now study the influence of the starting point of rays. We have the following result 

Proposition 3.1 Among oblique closed rays of angle a, the ray which spends the least time in uj 
in average is a ray which meets a point with integer coordinates . 
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Proof : Let B = {xo,yo) be a point of T such that the ray 7 — [B,a] does not meet a point 
with integer coordinates. Let us adopt the point of view of infinite plane and let Bq, Bi, B2, ■ ■ ■ 
be the intersection points of 7 with the lattice Z + Z. A direct consequence of Theorem 11.11 is 
that BoB2N{p+q) = L = 2N^Jp^ + q^. u ~ {pi + qj)/\/p'^ + is the direction of the ray. Let 
n = {—qi + pj) / \/ p'^ + q^ be such that (u, n) is a direct orthonormal frame. For a point P = {k, I) 
with integer coordinates, let us compute the algebraic distance to the ray 7 

d{P, 7) = ^.rt = P^-lk-{Py^-^^^) 

VP^ + q^ 

There exists a point A with integer coordinates and which verifies: 

7) = min{d(P, 7)|P has integer coordinates and d{P,j) > 0}. 




Figure 9: Case of a ray which does not meet any point with integer coordinates 

Let Ai be the intersection point of the line meeting A and directed by u (the points A solutions of 
the minimization problem above are clearly on a same line directed by u) and the line belonging 
to the lattice Z + Z which contains Bi (see Figure [9] for the case = 2). In the same way, there 
exists C with integer coordinates which satisfies 

c?(C, 7) = min{c?(P, 7)|P has integer coordinates and d{P,j) < 0}. 

Let us note Ci the intersection of the line meeting C and directed by u and the line belonging to 
the lattice Z + Z which contains Bi (see Figure [9]). 
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There arc three obhque rays of angle a 7^ = [AojO;], 75 — [B,a\ = [Bo,a\ and jc ~ [Co,Q^]- 
According to the definitions of A and C, the segments \Ai Aj+i], \Bi and [Cj Ci+i] belong to 

the same square Cfe,;. We define Sfc = 1 if [.Bfe -Bfe+i] C wt, fffe = if [B/j -Bfe+i] wt- Then, if 



m(7B) = y 



[-Bo-B2JV(p+g)] 



_ 1 



2N(p+q) 



we have: 



_^ 2Ar(p+g) ^ 27V(p+q) 

m(7A) = y ^ EkAk-iAk and m(7c) = j ^ ekCk-iCk- 




These three rays are parallel. Hence, there exists t e (0, 1) such that for all i, Bi is the barycenter 
of {Ai, t) and {Ci, 1 — t). If the quadrilateral Ai Ai+i Ci C^+i is a trapezium (or in the degenerate 
case a rectangular triangle), Thales' theorem implies that BiBi^i = tAiAi^i + {1 — t) CiCi+i. If 
the quadrilateral is a parallelogram, this inequality remains valid because the three lengths are 
equal. As a consequence, we have 

m(7B) = tm{'yA) + (1 - t)m{'yc)- 

The minimum of the three quantities 771(7^1), m{'-fB) and m{'-fc) is attained for one of the two rays 
J A ou 7c. This shows that 7 cannot be the minimum in the set of oblique closed rays. □ 

It remains to find the angles a wc need to consider among the rays which meet a point with integer 
coordinates. In the fiat torus T = [0, 2N]'^, there are 4iV^ possible starting points. However, the 
ray 7 = [{xo,yo),a] with {xo,yo) G meets 2 A'' points with integer coordinates. Indeed, 7 meets 
the point {xi,yi) if and only if there exists A € [0, 2N) such that 



Therefore, the points belonging to 7 and with integer coordinates are exactly the 2A'' points ob- 
tained for A = 0, 1, . . . ,2N — 1. There exist at least 2A' rays which meet a point with integer 
coordinates. 

Proposition 3.2 Let d be the CCD. of p and 2N , let d' he such that dd' = 2N . Then, the 2N 
rays of angle a starting at with < i < d— 1 and < j < d' — 1 meet one and only one time 
all the points of T with integer coordinates. 

Proof : Let p' be such that dp' = p. We note that p' and d' are relatively prime. Let also 71 
and 72 be the rays starting respectively at A and B with integer coordinates (fci, fe) and (fcs, k^) 
(with < fci, fcs < d — 1 and < k2, k4 < d' — 1 ). There exists a point C e 71 fl 72 with integer 
coordinates if and only if there exist four integer numbers A, ji, n and m such that 



Therefore, fci — fcs = (/x — X)p + 2Nn and fci — fcs is divisible by d. Since — (rf— 1) < fci — fcs < d—1, 

we have fci = k^ and (/x — X)p = —2Nn. As a consequence, (p — X)p' = —d'n and since p' and d' 
are relatively prime, (/it — A) = kd', where /c G Z. In the same way, k2 — k^ = {fi — X)q + 2Nm = 
kd'q + 2Nm. Hence ^2 — ^4 is divisible by d'. Since —{d' — 1) < k2 — < d' — 1, we have ^2 = k^ 



Xi = xo + Xp [2] 
yi=yo + Xq [2] 



X = ki + Xp = ks + ijp + 2Nn 
y = k2 + Xq = ki + fiq + 2Nm 



and 71 = 72. 



□ 
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3.4 The algorithm 

In this section, we give a method to compute expHcitly g when w is a finite union of squares Cij . 
At first, according to Propositions 11.11 and 13.11 the study can be restricted to closed rays starting 
at a point with integer coordinates. 

The user of the program must input the matrix M which represents w and the value of a parameter 
PQmax which avoids infinite loops (this problem never appeared until now). 

The algorithm is the following 

1. Compute the minimum g of the 2N numbers 

1 ^ 1 ^ 

-^M,,, and ^E^^.^' 

i=i j=i 

which corresponds to the minimum of m(^) when 7 is a vertical or horizontal ray. 

2. From the matrix Af , build the matrix A defined in 



3. For all (p, q) such that p>l,9>l,p + <z< 2A^ and p and q relatively prime, compute the 
number (the way to compute TOp,g is explained below) 



<'p,q 



mm 



{to(7)|7 = [X, arctang/p]}. (3.4) 



Xe{l,2,...,2Af} 

Then g <— minjg, mp^q}. 
4. If at this step g ~ \uj\, the user inputs the parameter PQmax which corresponds to the great 



est product pq which will be considered ; in the other cases PQmax = E{ . ''Y min{|cLi|, 1 — 



|}) + 1- 



\UJ 

5. Find all the (p, q) relatively prime such that p>l,q>l,p + q> 2N, pq < PQmax. Sort 
these couples from the lowest product pq to the greatest one. This gives a list L containing 
n couples. 

6. z <— 1 and as long as — 5 > min{|w|, 1 — |a;|} and i < n, 

(7 <— min{(7, mp.g} and + 
where {p, q) is the i-th couple of L and mp g is defined above (see ^.4^ ). 

7. Finally 

• if at the end of this loop, g = \lu\, then 



PQmax 



min{|t^|, 1 - \uj\} < g{uj) < \uj\. 



This means that one of the following cases occurs; g{u!) = \lu\, or too few families of 
rays have been considered. 

if g 7^ \uj\, according to Theorem 13.11 g{Lo) is equal to g. 
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We now present the method of computation defined in (|3.4p for two relatively prime integers 
P > 1, 9 > 1- At first, let us consider the ray 70 — [0, arctang/p] and let us study the way to 
compute 771(70). Its period is L = 2Ny^p'^ + and between the instants t = and t = L, it meets 
2A^(p + q — 1) + 1 points which have at least one integer coordinate. Let -P2Af(p+g-i)+i be 

these points (see Figure fTO]). Remember that ^';(p+g-i)+i, 1 < ^ < 2iV have integer coordinates. 



Pl7 





P 


p ^^y'^ 

p IS/^ 

14/ 

P13/ 


p 

1 







1 23456789 10 11 12 



Figure 10: Ray 70, for iV = 2, p = 3 and q = 2 
Let (ife), (jfc) and (sfc) be the three other finite sequences (1 < fc < 2iV(p + q — 1)) defined by 

PkPk+i C [ifc,ife + 1] X [jfc, jfc + 1] 

and 

^ f 1 if [ife, ife + 1] X [jfe, jfe + 1] C w 
\ if not 

Then, 777(70) = X]fc=i''"''' £kPkPk+i, the computation of 771(70) is equivalent to the computa- 
tion of the sequences (PkPk+i), (*fc) and (jk)- 

• Computation of the lengths PkPk+i 

At first, note that it suffices to compute the lengths PkPk+i for 1 < k < p + q — 1- Indeed, 
Pk+(p+q-i)Pk = pi + qj, hence the sequence [P^Pk+i) is periodic of period p + g — 1. Since p and 
q are relatively prime. Pi and Pp+q_i do not have integer coordinates. 

Then, the equations of the ray 70 are x ~ ptj \/p^ + q^ and y = qt j \J pp- ^ (p- . It meets the vertical 
lines defined by a; = 7 at instants — i\Jp^ + q^ jp at points [i,iq/p) and the horizontal lines 
defined by y = j at instants t'^ = 2\/p^ + (p I q and at points {jplq-,j)- 
Let us then define the list by 

- 1 < i < 9 - 1} U {U, 1 < 7 < p - 1}. 
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Let be the list obtained from by sorting its elements in increasing order. Let also a be the 
permutation of {1, 2, . . . ,p + q — 2} such that 

Let and be defined by 

= {0} UV'U {V^T^} and = {Vi - V,^ Vi -Vi,..., V^+^ - J. 
Finally, let V be defined by 

V= V*[JV*U---[JV* 



2N times 

then for all fc, 1 < fc < 2N{p + q - 1) PfePfe+i = Vk and 

2N(p+q-l) 



2NsJv^ + q^ fr[ 



As an important remark, it may be convenient to work with pq/ y'p^ + q"^ .V instead of V. Thus, 
^1 = {pjj 1 < j < 9 - 1} U {qi: 1 < * < P - 1}, ^3 = {0} UV^ U {pq} and hence, the elements of V 
are all integer numbers. m(7o) is then given by 

= 2Npq • 

As a consequence, m{jo) is a rational number whose denominator and numerator are explicitly 
known. 

• Computation of the sequences (ik) and (jk) 
Let and be the two following lists 

/!= {1,1,. ..,1, 0,0,. ..,0} andJi= {0,0,. ..,0, 1,1,.. .,1} 

^ V ' ^ V ' ^ V ' ^ V ' 

q — 1 times p — 1 times (Z — 1 times p — 1 times 

This means that the numbers 1 which appear in the sequence correspond to the points Pk 
which belong to horizontal lines and the numbers 1 of the sequence correspond to the points 
Pk belonging to vertical lines. From these two lists and from the permutation a, we define 

= (^{^^) = {-^ct(1)'-^ct(2)! • • • i-^CT(p+g-2)} and = Cr(J^) = {-/^^(i) , -^^^(2) , • • • , Ja{p+q-2)} 

and also 

[- 2 P+q-2 ^ f 2 P+q-2 ^ 

0,/2,^/2,..., Il\ and J3= 0,J2,^J2,..., 

I k=l fe=l ) I fc=l k=l ) 

Finally, we set 

/ = /3 U (p + J3) U • • • U {[2N - l]p + P) and J = J3 U (g + J3) U • • • U {[2N - l]q) + . 
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Then, for all fc, 1 < fc < 2N{p + q ~ I), iu = I + h and j ^ = 1 + Jfc. 

Let us introduce the matrix A defined from M in (|3.3p and a function r defined over integers by 



r{k) G {1, 2, . . . , 2N] and 3l^1\k = 2N.I + r{k). 
Then, for all fc, 1 < fc < 2N{p + q - 1), Sk = Aiik).rUk) = A(i+ik).r(i+Jk)- 

There are two reasons why we consider all rays of angle arctang/p at the same time in the compu- 
tation of mp g. First, all these rays have the same sequence Vk- In addition, let 7 — [Xq, arctan(^)] 
where Xq = (a;o,yo) has integer coordinates. Then, the sequences (ik) and (jk) of 7 can be com- 
puted from the lists Ik and Jk by the formula: ik = 1 + xq + Ik and jk = I + Xq + Jk- 

For two relatively prime integers p > 1 and q > I, the number rupq can be computed in the 
following way 



1. compute the CCD. d of p and 2N and from Proposition 13.21 the 2N starting points with 
integer coordinates that must be considered are known. Let us note Xi,X2, ■ ■ ■ ,X2n and 
Yi,Y2, ■ ■ ■ , Y2N their coordinates. 

2. Compute the three lists of 27V (p -t- g — 1) elements Vk, Ik and Jk defined above. 

3. TOp.g is then the minimum of the 2N numbers 

2Npq 

All these quotients have the same denominator and hence, are easy to compare by looking 
at their numerator. Finally, this gives the explicit value of g. 



We now present the results obtained with our algorithm for explicit domains uj. 



area = 1/2, 401 rays, g = 1/3 




40 60 
product p*q 



Figure 11: An example of computation of g{Lu) 
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Figure 12: Minimizing rays : m{'y) = g{uj) = 1/3 




As an example, for the computation of g{u}) for the domain oj of Figure [71 the minimum for 
vertical and horizontal rays is g = 1/3 and there are 45 families of oblique rays of first kind (i.e. the 
rays of part 3 in algorithm). Among these rays, there exists 7 such that rn{'y) = 1/3 but there is no 
7 such that m(7) < 1/3. At this step, the value of g is 1/3. There are also 344 famihes of obHque 
rays of second kind (corresponding to part 5 of algorithm). None of them satisfies 771(7) = V^- 
We obtain g{uj) = 1/3. 

Figure [11] presents the values of mp^q as a function of pq and presents also inequality of Theorem 
13.11 Note that this inequality is not optimal for this particular domain. However, (|3.2p is optimal 
for a unique square Cij. This shows that the error comes from the summation of inequalities (j3.2p 
for all the squares ofujT- Nevertheless, it is useful because it certifies that the value of g{uj) is exact. 

There exist 10 rays which satisfy m{'j) = g{uj); six vertical or horizontal rays (more precisely six 
families of vertical or horizontal rays). They are represented on Figure [T2la. There are also two 
rays of angle 7r/4 represented on Figure [T2lb, one ray of angle arctanl/2 represented on Figure 
[T^c. and also the ray 7 = [0, arctan 1/4] represented on Figure [T^d. 

The domain oj of Figure [T3l of area |a;| — 29/100 is represented here to show that intuition can 
be false. Indeed, consider only vertical rays, horizontal rays and closed rays of small period. The 
intuition says that g{uj) = 1/5. However, the ray 7 = [0, arctan8/7] is such that 771(7) = 111/560 
(note that 1/5=112/560), that is the exact value of g{uj) as shown by the complete computation. 

4 Appendix : proof of Theorem 12.1! 

First, it is clear that 

5(^)<5'(^)- (4.1) 
Indeed, for all 7' G C and all t > 0, m{'y^) > inf-ygc 'm{'yt)- Hence, 

777(7') = limsup777(7j) > limsup inf m^jt) = gi^)- 

t-»+oo 76C 



Since this equality is true for all 7' G C, this shows (|4.ip . 
Let us prove that 

g{Lu) > g\uj). (4.2) 

Let us fix e > and choose a sequence (7")n C C and a sequence {tn)n which tends to +00 such 
that 



hm r77(7rj=5(^). (4.3) 

For all 77 e N, there exist X„ G and 0„ G [0, 27r[ such that 7" = [Ar„, 0„]. Up to a subsequence, 
we can assume that there exists {Xq, 6*) G x [0, 2tt] such that 

VnnXn = Xq and lim0„ = 9. 
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We claim that there exists L > such that 

m{[v,e]L) > g'{u) - e for all ye^. (4.4) 
If tan(0) G K\Q, then set L — to, where to is as in Proposition ll.il Indeed, together with relation 

mi[y,e]L)> \uj\-e>g'{uj)-e. 

If tan(0) e Q or if ta,n{6) — oo, the ray [y, 9] is periodic and the period T does not depend on 
y Hence, 

m{[y,9]L)^m{[y,e])>g'icj). 

Thus, let us set L ~ T. Clearly, t„ can be assumed to be a multiple of L. Indeed, replacing f„ 
by tn + Sn (where |s„| < L) does not change the limit of m{j" ). Now, assume that 0„ = 9. If 
tan(a) G Q or if tan(Q;) = oo, then clearly lim„ m(7" ) — m(j") = 0. If tan(Q;) G this relation 

follows from Proposition ll.il Then, by relation (|4.3|) . it follows that 

g{uj) > liminfm(7") > g'{uj). 

n 

This proves (|4.2p . Hence, we can assume that 

9n + 0. (4.5) 

Now, let us define, for k G {0, ^-1}, Xk,„ = 7"(/cL) and -/>''"{t) = -/"(kL + t) for t G [0,L]. Let 
us define also 7'='"''* as the unique trajectory of length L and slope 9 such that ^'=."■.''(5) = 7'='"(s) 
and 

se[o,L] 

. k.n.N 
7 

7 ■ ■ 7 (s) = 7 (s) 

^fe,n,0 

Now, we prove that 

Step 1 There exists a sequence {kn)n such that 

|^fc„,n n 9w| < e and 771(7''"^") < g{w) + e 
for all n large enough. Here, |^fc„,n H dwl denotes the length of the curve ^fc„,n H dw. 
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In the following, the length of any curve C will be denoted by |C|. Let 



/C„ = <j A; e {0, ^ - 1} s.t. m(7'^^") < g{w) + e 



Set 



and 



M= [\\ U ^^^^ 

mGA \ n>m ;A;GACti 



n U 

For X G A^, let 7x denote the trajectory of length L and slope 9 such that 7x(0) = X. Note that 
and M are not empty. Otherwise, /C„ would be empty too and we would have 771(7*^'") > g{uo)+e 
for all n > uq and forall fc. Thus, we would have 771(7" ) > 9{^) + £ that contradicts (|4.3p . Now, 
we prove that 



- U 7,- (4.6) 

Take a sequence (X/^_„)„ such that lim„X;^,„ = y G A4. Then, for all t G [0,L], lim7'"'"(t) = 
7j,(i). We obtain that 

U c ^- 

Conversely, if z G then there exist two sequences (/«)« such that Z„ G ICn and (s„)„ C [0, L] 
such that z = lim„ 7'"'"(s„). Up to a subsequence, we can assume that there exist y ^ Ai 
and s G [0, L] such that lim„Xfe^ .„ = y and lim„ s„ = s. As one can check, this implies that 
z = lim„ 7'='»'"(s) = 7y(s) and therefore 

yeM 

This proves (|4.6p . 

We now distinguish two cases. 

First, assume that there exists an infinite sequence {yn)n in M. such that 7y^ '~'7iy„/ = if f^ 7^ 
Since |(9a;| is finite, one can choose y G (7/ is one of the 7/„) such that \^y n (9a;| < |. Otherwise, 
for all 71, \duj fl 7y„| > | and hence, since 7y„ fl 7^^, = 0, we would have 

+oo = ^|atjn7yJ = |(u„7yjnacj| < \du\. 

n 

This is false. The definition of M implies that there exists a sequence (fc„)„ such that y = 
lim„Xfe^^„ with fc„ G /C„. Clearly, for all < G [0,i], 



lim7'="'"(t)=7y(i). 
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(4.7) 



Set Bn = Um>ri ^k„,n- Obviously, Bn is a decreasing sequence of sets. We recall that 

se[o,L] 

By construction, for all t e [0, L], = 7y(i)- This imphes that HneN — 7a- Since 

l7j, n Swl < |, for all n large enough, we have |i3„ n duj\ < e. Since ^fe„.„ C i?„ and since e /C„, 
this proves Step 1 in this case. 

Now, assume, for all infinite sequences {yn)n in there exists n ^ n' such that 7y„ ^ly^, 0- 
Clearly, this case occurs if M is finite. Then, choose a finite set • • • , j/fe} C M. such that 
ly-n ^ ly„' = if ?^ 7^ and assume that k is maximal. This means that for all ?/ G A^, there 
exists rt S {l,...,fc} such that 7j, fl 7y„ 7^ (otherwise the set {y,yi,--- ,yk] would have the 
above property and k would not be maximal). Since |7y| = L, we have: 7j, C 7^^ where 7^^ is 
the trajectory of length 3L + 2 and slope 9 such that l'y^{L + 1) = 7j/„(0). Note that instead of 
7^^, one could have considered trajectories of length 3L. We need to consider trajectories strictly 
greater than 3L to get relation (|4.10p below. For simplicity, we choose 3L + 2. Using (|4.6p . we see 
that 



^^[jlL- (4-8) 



n=l 

Let n be fixed. Let 



7"(t) G U 



m— 1 



The set Um=i 7ym ^ finite union of disjoint trajectories of finite length and slope 0. Moreover, 
^ji 7^ ^ (see (j4.5p ). Hence, A„ is a discrete set of points. Since Um=i 7i/„, is closed, is closed 
too. This shows that A„ is finite. Let now ti,t2 G ^„ with ti < t2. As one can check, there exists 
a > independent of n such that if t < a\9n — ^|~^, 7"(*i + ^) ^ Um=i 7i/™- Indeed, lJm=i 7y„ 
is a finite union of segments of slope 0. Then, it is clear that the trajectory 7„, whose slope is 
meets two different segments of lJm=i ly^ ™ ^ time proportional to the difference |0„ — 0| of the 
slopes. Hence t2 — ii > a\6n ~ 0\ ^. This proves that 

#An<-\6n-0\. (4.9) 
a 

Moreover, 

m(7,:) = - / x^{ii{t))dt = J2t WMdt - r E '"(7'")- 

From the definition of Kn, we write that for I G /C„ we have 771(7''") ^ ^ind we know that for 
I e {0, - 1}\/C„, we have 771(7''") > ^(w) + e. We obtain that 

m(7rj > ^# ({0, I - 1}\/C„') (g(c.) + e) = ^ - #/C„') (g(c.) + e). 
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Since lim„ m(7j^) — g{uj), t 



where c > is independent of n. It follows from (|4.9p that we can choose /c„ e /C„ such that 
[fc„L, (fc„ + l)L]nA„ = 0. Hence, 7''"^"nljfjj^;^ 7^,^ =0. We can assume that lim„ Xfe„,„ = ?; G A^. 
Then, for all t G [0,i], we have: lim„ 7'^"'"(<) = [y,d]{t). Since 'y'^^ has been chosen such that 
their length is 3L + 2, there exists an m such that 

[y,0]LCj'y^ (4.10) 

with [y, 0](O) 7y„, (0) [y,9]{L) ^ jy^^{L). As one can check, this implies that ^fc„.n H 
Um=i 7j(„ = ^ (^66 figure below). 



yen 



7" 



7. ^'^^^ 



i 

lym 

By (|4.7p . we obtain that 



Now, set 



and 



Ak^^nnM=9. (4.11) 
dn=dist{Ak„^„,j' )= inf dist{X,j' ) 

= sup dist{X,j' ). 

We have < d„ < rf^. Since the sequence of sets {Ak,^,n)n converges to 7y (see figure above) and 
since "fy C 7y^, it is clear that lim„ d'n = lim„ (i„ — 0. After passing to a subsequence, we can 
assume that for all n, rfjj+i < dn- In other words, the sets Afc„,n can be assumed to be all disjoint. 
Since \du!\ is finite, there exists a infinite number of ^fc„,n such that n dLu\ < e. We can 

assume that, for all n, \Ak^^n H du}\ < e. Since /c„ e /C„, this proves Step 1 in this case. Note that 
the sequence (fc„)„ does not necessarily tend to +oo. 

Let now 7" — 7*-'"'". Keeping the same notation as in the definition of Ak^n- Furthermore, let 
7o = -y^-'^.'^.o. This means that 7g is the unique trajectory of slope 6 and length L such that 
7^(0) — 7"(0). We now prove that 
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step 2 We have 

M7")-rn(7o")| <|e 

for any n large enough. 

The rays 7" and 7g can be seen as segments of length L in the plane. Let p„ be the projection 
on 7q in the direction of the line (^"■{L) 7q (L)). Thalcs' theorem implies that for all t e [0,i], 

Pn(7"(t))=7o"W- 



7oW=Pn(7"(i)) 
7"(0)=7o"(0) 

Now let t e [0, L]. Assume that A'„(7j(t)) ^ X^{jQ{t)). Then, there exists a point X of duj between 
7"(0 and7ff(i). Thus, X e [7"(t), 7j(t)] f] Sw. This shows that p„(X) G p„(5w f] Afe„,n)- Indeed, 
remark that, for all t G [0, L], 

[7"W,7o"W]c^fe„,„. 

Since 70 is parametrized by the length, the segments [0, L] and 7q can be identified. We obtain 
that 

{t G [0,L] s.t. ^-^(7"^) 7^ ^-^(70 W)} Cp„(ac^f|^fe-")- (4.12) 

Now, write that 
By (|4T2)) . this gives 

l™(7")-"i(7o)l < ^|p„(ac.f|Afe„,„)|. 

Remember that p„ is a projection whose angle tends to ^. Hence, for any piece of curve C, 
bri(C)| < 2|C| (the constant 2 could be replaced here by a constant c„ which goes to 1 with n). 
It follows that 

|m(7")-m(7o")| <||ac.f|Afc„,„|. 
Step 2 is then a direct consequence of Step 1. 
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step 3 Conclusion 

The trajectory 7g has length L and slope 9. The definition of L (see (|4.4I) ') then implies that 

™(7o ) > 5'(^) - 

By Step 1, 

m(7") < .gH + e. 

Step 2 then shows that 

2 

g'{uj) - e < g(t^) + e+ -e. 
Since e is arbitrary, we obtain (|4.2p . Together with (|4.ip . the proof of the theorem is complete. 
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